home *** CD-ROM | disk | FTP | other *** search
/ IRIX Base Documentation 2002 November / SGI IRIX Base Documentation 2002 November.iso / usr / share / catman / p_man / cat3 / SCSL / dlarrf.z / dlarrf
Encoding:
Text File  |  2002-10-03  |  4.6 KB  |  133 lines
  1.  
  2.  
  3.  
  4. DDDDLLLLAAAARRRRRRRRFFFF((((3333SSSS))))                                                          DDDDLLLLAAAARRRRRRRRFFFF((((3333SSSS))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      DLARRF - the initial representation L D L^T and its cluster of close
  10.      eigenvalues (in a relative measure), W( IFIRST ), W( IFIRST+1 ), ..
  11.  
  12. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  13.      SUBROUTINE DLARRF( N, D, L, LD, LLD, IFIRST, ILAST, W, DPLUS, LPLUS,
  14.                         WORK, IWORK, INFO )
  15.  
  16.          INTEGER        IFIRST, ILAST, INFO, N
  17.  
  18.          INTEGER        IWORK( * )
  19.  
  20.          DOUBLE         PRECISION D( * ), DPLUS( * ), L( * ), LD( * ), LLD( *
  21.                         ), LPLUS( * ), W( * ), WORK( * )
  22.  
  23. IIIIMMMMPPPPLLLLEEEEMMMMEEEENNNNTTTTAAAATTTTIIIIOOOONNNN
  24.      These routines are part of the SCSL Scientific Library and can be loaded
  25.      using either the -lscs or the -lscs_mp option.  The -lscs_mp option
  26.      directs the linker to use the multi-processor version of the library.
  27.  
  28.      When linking to SCSL with -lscs or -lscs_mp, the default integer size is
  29.      4 bytes (32 bits). Another version of SCSL is available in which integers
  30.      are 8 bytes (64 bits).  This version allows the user access to larger
  31.      memory sizes and helps when porting legacy Cray codes.  It can be loaded
  32.      by using the -lscs_i8 option or the -lscs_i8_mp option. A program may use
  33.      only one of the two versions; 4-byte integer and 8-byte integer library
  34.      calls cannot be mixed.
  35.  
  36. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  37.      Given the initial representation L D L^T and its cluster of close
  38.      eigenvalues (in a relative measure), W( IFIRST ), W( IFIRST+1 ), ... W(
  39.      ILAST ), DLARRF finds a new relatively robust representation L D L^T -
  40.      SIGMA I = L(+) D(+) L(+)^T such that at least one of the eigenvalues of
  41.      L(+) D(+) L(+)^T is relatively isolated.
  42.  
  43.  
  44. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  45.      N       (input) INTEGER
  46.              The order of the matrix.
  47.  
  48.      D       (input) DOUBLE PRECISION array, dimension (N)
  49.              The n diagonal elements of the diagonal matrix D.
  50.  
  51.      L       (input) DOUBLE PRECISION array, dimension (N-1)
  52.              The (n-1) subdiagonal elements of the unit bidiagonal matrix L.
  53.  
  54.      LD      (input) DOUBLE PRECISION array, dimension (N-1)
  55.              The n-1 elements L(i)*D(i).
  56.  
  57.      LLD     (input) DOUBLE PRECISION array, dimension (N-1)
  58.              The n-1 elements L(i)*L(i)*D(i).
  59.  
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. DDDDLLLLAAAARRRRRRRRFFFF((((3333SSSS))))                                                          DDDDLLLLAAAARRRRRRRRFFFF((((3333SSSS))))
  71.  
  72.  
  73.  
  74.      IFIRST  (input) INTEGER
  75.              The index of the first eigenvalue in the cluster.
  76.  
  77.      ILAST   (input) INTEGER
  78.              The index of the last eigenvalue in the cluster.
  79.  
  80.      W       (input/output) DOUBLE PRECISION array, dimension (N)
  81.              On input, the eigenvalues of L D L^T in ascending order.  W(
  82.              IFIRST ) through W( ILAST ) form the cluster of relatively close
  83.              eigenalues.  On output, W( IFIRST ) thru' W( ILAST ) are
  84.              estimates of the corresponding eigenvalues of L(+) D(+) L(+)^T.
  85.  
  86.      SIGMA   (input) DOUBLE PRECISION
  87.              The shift used to form L(+) D(+) L(+)^T.
  88.  
  89.      DPLUS   (output) DOUBLE PRECISION array, dimension (N)
  90.              The n diagonal elements of the diagonal matrix D(+).
  91.  
  92.      LPLUS   (output) DOUBLE PRECISION array, dimension (N)
  93.              The first (n-1) elements of LPLUS contain the subdiagonal
  94.              elements of the unit bidiagonal matrix L(+). LPLUS( N ) is set to
  95.              SIGMA.
  96.  
  97.      WORK    (input) DOUBLE PRECISION array, dimension (???)
  98.              Workspace.
  99.  
  100.      IWORK   (input) INTEGER array, dimension (???)
  101.              Workspace.
  102.  
  103.      INFO    (output) INTEGER
  104.              = 0:  successful exit
  105.              > 0:  if INFO = 1, IFIRST and ILAST had illegal values
  106.  
  107. FFFFUUUURRRRTTTTHHHHEEEERRRR DDDDEEEETTTTAAAAIIIILLLLSSSS
  108.      Based on contributions by
  109.         Inderjit Dhillon, IBM Almaden, USA
  110.         Osni Marques, LBNL/NERSC, USA
  111.  
  112.  
  113. SSSSEEEEEEEE AAAALLLLSSSSOOOO
  114.      INTRO_LAPACK(3S), INTRO_SCSL(3S)
  115.  
  116.      This man page is available only online.
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.